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163. 


ON HANSEN'S LUNAR THEORY, 


[From the Quarterly Mathematical Journal, vol. 1. (1857), pp. 112—1235.] 


THE following paper was written in order to exhibit, in as clear a form as may 
be, the investigation of the remarkable equations for the motion of the moon established 
in Hansen’s “Fundamenta Nova Investigationis Orbite vere quam Luna perlustrat," 
&e., Gothe, 1838. I have availed myself for this purpose of the remarks in Jacobi's 
two letters in answer to a letter of Hansen's, Crelle, t. xri. [1851], p. 12; it may 
be convenient to remark that the quantity there represented by A, and which does 
not occur in Hansen's own investigation, is in this paper represented by 6. 

The position of the moon referred to the earth as centre is determined by 


r, the radius vector, 
L, the longitude, 
A, the latitude. 
Suppose, moreover, that the attractive force at distance unity, — x (M + E), is represented 
2548 
by na? then the principal function will be y--7, and the disturbing function R 


may be represented by n’a*Q; the expression for the half of the vis viva is 


rfe («e (0) 


and the equations of motion are therefore 


d dr T S? A (Sy 3) n'a? = nma’ dQ 
ae ea dt nna po ar’ 
3 (r cos? A A = na? 7 s 
d dA : dL? 4.5 04 
q (rg) e eos A sin A (7) = na T E 


where Q is considered as a function of r, L, A. 


www.rcin.org.pl 


14 ON HANSEN’S LUNAR THEORY. [163 


Consider the orbit as an ellipse, then putting 
a, the mean distance, 
n, the mean motion = "| CE , 
e, the excentricity, 
c, the mean anomaly at epoch, 
w, the distance of perigee from node, 
0, the longitude of node, 
i, the inclination, 
c, the distance from node, 
W, the reduced distance from node, = L — 0, 
U, the excentrie anomaly, 
f, the true anomaly, = e — o, 


the elements of the orbit are a, e, c, e, 0, à, and we have from the theory of 
elliptic motion, 
ntt+c= U-—esin U, 


V(1 — e) sin U 


ison - 
dmn cos U—e ' 
a (1— e) 
«a (l+ SALA C ey 
r=a(1—ecos U) VIP 


Moreover i is the angle at the base of a right-angled spherical triangle, the base, 
perpendicular, and hypothenuse of which are WV, A, P, hence- 


tan V = cos? tan  , 
sin A —sinzsin  , 
sin Y —cot;tan A , 
cos ® = cos A cos Y. 
Considering the elements as constant, we have 


E _ maesin f 


~ Va =e)’ 
K, na? y (1 — e) 
T? 
2 na? /(1 — e) 
dt ri : 
dv _ cost na? (1 — se 
dt cos? A g 


dL _ cost na’ y(l -— e) 
dt cos? A 7? 


na? /(1 — e) 
T? j 


on sin 4 cos V 
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Hence also 
d (5) n?a*e cos f 
dt \dt TC CERES c 
d dL 
di (e cos? A 23 = 0, 
io) see ee) 
dt\ dt cos? A ri > 
d (df Ma.. 
AFA =- grean j 


The foregoing values show that the equations of motion, neglecting the terms which 
involve the disturbing functions, are satisfied by the elliptic values of r, L, A: and 
in order to satisfy the actual equations of motion, we have only to consider the 
elements as variable and to write 


dr =o, 

dL = 0, 

dA. = 0, 

d = na? = dt, 

d (v cos? A T) = ma? > dt, 
dA dQ 

d (r a) == ng? dA dt, 


where the differentiations relate only to the elements, or, what is the same thing, to t 
in so far only as it enters through the variable elements: the system is at once 
transformed into 


dr zz, 
dL = 0, 
dA =0; 
nae sin f oe UU 
d — &) = niat 3. dt, 
d na? J(1 — &) cos = nta ae dt, 
dQ, 


d na? /(1 — &) sin cos V = na? EX: dt. 
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Now W-L-—6, or (supposing, as before, that the differentiations relate to f, only in 

tan ® 
di; 

sin 4 


so far as it enters through the variable elements) dV = — dó, and thence d0 = 


we have also db = — cos? d0. The equations containing e and e give 


cos 4 d na? /(1 — e) — na? y (1 — e) sin i di zu zn dt, 


cos V sin 4 d na? /(1 — e) + na? /(1 — e) (cos V cos ¢ di + sin V sin 4 d0) = wa? e dt ; 


or, expressing dô by means of di and reducing, the second of these equations becomes 


ee di poi d dt, 


* LI 1 
sin id na? J(1 — e) + na? (1 — €) — cosy "co V dA 


and combining this with the first of the two equations, and observing that 


cos? ù a^ TES 1 
cost A cogi sp * SI * = eor: 
we find 
dna?4(1—28)-— na’ fum = + sin 4 cos V A) dt, 
1 na? dQ dQ, 
di irm pao de sin 4 cos? Y 27. + cos i cos Y 5) dt. 


Now, considering Q as a function of r, 0, 7, b, then A, L are given as functions of 
0, i, ® by the equations sin A=sinisin®, tan V = cositan b, Y —L—60, and after 
some simple reductions, 


do dQ 
dr dr’ 
do dg 
dó aL’ 
dQ dQ dQ 
q= m e(- sin 4 cos? y + cosi cos Y a). 
dQ, _ cosi dQ dQ 
m" ( SUA dr nice YF) 
whence also 
dQ oc — in 4 rom 
dq) d ~ Sin tco di' 
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We have therefore 


d na? (1 — &)= nas dt, 
di A ram S IT 
de = ERE = dt, 

do - vale qi 

dr EE 0, 

d Tm a = nas > dt. 


Suppose now that we have 
p, a radius vector, t for t, 
$, a true anomaly, do., | 
v, an excentric anomaly, do., 
Le. let p, $, v, be what the radius vector, the true anomaly and the excentric 


anomaly become when the time £, in so far as it enters directly, and not through the 
variable elements, is replaced by a new variable 7. We have 


"T --c—v-—eSsinv, 


$ ulia (1 — &)sinv 
cosv—e ’ 

y . a(1—e). 

p — (1 f conu) = Ta Cond" 


and of course the differential coefficients of p, $ with respect to r may be at once 
deduced from the corresponding expressions for the differential coefficients of r, f with 
respect to f£, the elements being considered as constant. Now, using l to denote a 
logarithm, and supposing that the differentiations affect only the elements, we have 


dla = da  2ede _pcosġde , pesin ddd 

P PET ET T T ' 
da  2ede _rsinfde , resin fdf. 
a l-e a(1—-e) a(1—e) 
and putting for shortness 


Simul pesin ddd 


a(l- e)’ 
_ 7, Te sin fdf 
X =dlr vert 
we find 
.psinó. (|, psing\/da_ 2ede\ ^ p _ sin $ cos f 
: rain X (1 yd 1-é) a(1—e) (cos ¢ sin f ) de, 
C. III. 3 
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or reducing 


p sin $ $ y 
iae r sin RIT 
da _ Rus 


l \(sing- sin $ + esin(f— ¢)) ( la ra) REL da}. 


sin f (1 + e cos $) 
Write for a moment 


P=neV(1-@&), Q= nae wn R=% oer 


|" X(1-e)' 1-4 ecos f" 
so that 23 
a (1 ges e) = Tig? , 
Se ee 
e 3 P+ ur Rm nu R 
We have therefore 
da  2ede 2dP 2 
«i-e P "aü-e) 
1 2 um | 4M 2o BE 
ign (aa Pes R wa? x) ebd i; PQUQ* fa mai T wat B) ah, 
which, after reduction, becomes 
E 1 " x Ape _ (1r ecos f cos f 
de = nd é) (e sin? f+ 2 (cos f +e cost f)) d P4 x 4 (1 — e) sin f dQ MUTO IU ve dk, 


and substituting these values, 


X,- CES X= aang (2-200 (f- 4) + esin sin (f 4) l5 dP 


—a(1—e)sin(f-— ¢) waa" 


4, Gtecos f) cot f sin (f=¢) 1 
T Sacer. 0 - Id aR} , 


or substituting for X, X,, P, Q, R, their values 


| psinó pe sin $ Ss 
pio -— nae gj - diu 


zu —e) (2- 2 cos ( f — $) + e sin fsin(f — 2) 2 I$ 8j dna? /(1 — e) 


ithe o 1c nae sin f 
p sin (f V ces e? Tae) 
l- 

allem 99 OE a Varma er 
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Now 

die (2-2 008 (F— $) +e sin fsin (f €); 
therefore 


p sin $ pe sin $ "4 


1 
- eos (f-$)-1+ zat (cos f- 9) -1) cer s d a - e) 


1 nae sin f 
-psin(f-4) d —44J-é) 


+ PENCIL) cot fein (— 9) p wy dr. 


So far the variations of the elements have, in fact, been treated as independent ; 
naesin f . 


but if we substitute for dna? (1 — e), d Ve)’ dr, their values in the disturbed 
motion, the equation becomes 
pe sin $ = MERC men p M na dQ 
dlpt Ta —e) (9f d$) = (eos C£ $) L+ qo eU $) Ja ig 
na  dQ 


Consider now the point in which the orbit is intersected by any orthogonal 
trajectory to the successive positions of the orbit, or to fix the ideas, the orthogonal 
trajectory passing through T, the point in question may, for want of a recognised 
name, be called the “departure point;" and the angular distances in the orbit measured 
from this point may be termed “departures ;” the expression, “the departure," is to be 
understood as meaning the departure of the moon. Write now 

. X» the departure of the perigee, 
v, the departure, =/+ x, 
v, the departure of the node, —X-— ow, 
©, the longitude in orbit of departure point, = 0 — c. 

It should be remarked that y is not properly an element, ie. it is not a function 
of a, e, c, œw, 0, i without t, and in like manner o and 6 (which depend upon y) are 
not elements. 

We have from the geometrical definition 


dy — do + cos i d6 ; 
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and therefore 
do = cos i dé, 


d = (1 — cos?) dé. 


Moreover v, =® +ø, which gives (assuming that the differentiations are performed with 
respect to f£, only in so far as it enters through the variable elements) dv, — db + de 
=d®+cosidé, i.e. dv,=0, an equation which might have been assumed for the purpose 
of defining the departure point; the equation, in fact, expresses that the departure 
v, is measured from a point not actually fixed, but such that the increment of v, in 
the interval of time dt is the angular distance between two consecutive positions of the 
moon. : 


We have, as above noticed, de —cosid6, and thence and from what has preceded 


PS E 
RUTE X. 
TE na coti dQ 

~ Jüa-éà) di’ 


Now the position of the moon can be determined by means of the quantities 
r, v, OG, c, i; hence Q (which has been considered as a function of r, ®, 0, i) may, 
if we please, be considered as a function of r, v, ©, c, ? and from the differential 
relations 


dr =dr, 

dv, =d® + cos i dé, 
dO — (1 — cos 2) dé, 
do = do + cos i d0, 


di — di, 
we find 

4o do 

dr dr’ 

dO _ dQ 

db dy,’ 

dn 

qg esie) DE 
dQ dQ 

di di 


we have therefore 
dQ 1-— cos; dO dO 1 dQ 


We? Gee ae Tab werd 


or in virtue of a preceding equation 


dQ 1—cosi dQ 1 dQ 
p ek go - tanicot d dij 
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and effecting the substitutions, and collecting the results, 


d na? V/(1—e*) = mna? = dt, 


dv, 
dr = 0, 
d ween = war ay, 
eee a 


where © is considered as a function of r, v, 8, c, i 


Instead of c, i we may introduce the new quantities p, q defined by the equations 
p=sin ising, 
q =sinicosa; 

this gives sin?¢=p*+ q’, c — tan3P and retaining in the formule the sine and cosine 


of i, to avoid the introduction of irrational functions of p*?+q*, we have 


dO — (1 — cos?) dé = — — e t de = [og nd .(qdp — pdq), 


cos $ cos 7 sin? 4 


... gqdp — pdg 
iae see FL a: +cos 1%)’ 


which determines © by means of p and g. We have moreover 
dp =  sinZcosc doe + cos isin c di, 
dq =—sinisin ø de + cos $ cos o di, 


dQ _ sing dQ , cosa dO 
dp cosi di sini do’ 


dQ cose dQ sinc dQ 


dq cost di sin? da’ 


from which equations and the foregoing values of di and do we find the values of 
dp and dq; the other equations of the system remain unaltered, and we have 
therefore 


d na? J(1— e) = na? ES dt, 
dr =| Q, 
naesinf —— dO 
Es vesc nas d; dt, 
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_ macogi (dQ _ p. , dQ 
NS v-e) es ~ cosi (1 + cosi) 26) i 
|. ma cos'i do q dQ 
as ~ V1 =e) E * cosi (1 + cos 1) Ta) id 


where Q is considered as a function of r, v, 8, p, q. The symbols T 3 T as 


employed by Hansen, mean that the differentiations are to be performed as if © was 
a function of p, q, such that 


18 iy 4 2È dg- WPi . 
2 dp = + ag 1 cos i (1 + cos 2)’ 


the last two equations being therefore nothing else than what Hansen represents by 


na cos dQ 


dp = Je) dg ^^ 
dQ 

d na cos? i dt. 

270-6) dp 


Write now 
à, the departure, 7 for t, 


ie. X is what v, becomes when f, in so far as it enters explicitly, and not through 
the variable elements, is replaced by the new variable +; so that, in fact X—-Fx. 
The values of r, v, could be at once found from those of p, X by changing 7 into t; 
and to determine the values of p, A, Hansen proceeds as follows: 
writing 

^ —H (G5 t), 


lp =T (5, t) - B, 


where ¢ and 8 are new variables functions of T and ż, and II, T are arbitrary 
functional symbols; so that if z, w are what £, Æ become when r is changed into t, 
we should have 


v, — II (z, t), 
-I(t-H4w; 
then the foregoing equations give 
d ; dt 
= oF 
pn d 
-T'(& 055 eT, t, 
x m dB 
dr =T" (& t) 3- ur , 
di ; 
de TY (f T (9445; 


dt ’ 


[where the accents and strokes denote differentiation in regard to £, t respectively]. 
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Hence eliminating II’(g, t) and I"(£, t) we have 


dx d£ dxdt dt 
jd ad eu 


dipdg dipdt dB dt dB dg 


dg 
dud ae ok M Oe 


or, what is the same thing, 


dg dX 
dt dt 1 
dr dt dr 
dt dA dlp 
dB dB dt dip dt dip dr 
eee ale + II, (5 on r,(£, t), 
dr dr dr 
or writing 
dg dg 
pat dt Q.4d8 dg dt 
“dr de’ " d dr de’ 
dr dr 
we have 
d?X d?X d£ 
I dép. cde Gba: dr? dr 
T= dk aide 7d di * I, 6 Orga LE OR, 
dr (2) F3 dr 
dlp dip 
_dlp drd dr _ 
A PEE om r, (g t) 
dr dr 


Now from the equation »=¢+x where x is independent of r, the differential 
coefficients of X and p with respect to 7 are at once deduced from those of f, r 
with respect to t, and we have, 


oF tear e sin 

(d UN $, 
dà ma /(1 — e*) 
dr /— p? 4 
dlp nae sin $ 


dr. pvü-ey 
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Consequently | 
d mna?4(l—6&)  2pesin d dX 
~ na?^4(1—&) dt p? a(l— e) dt 


. 2pesin $ 
"a(1—e) wow na? es e) 


2 2pe sin $ dX 1 d 
a(l- e) dt t no? V(1 — e) dt neo 


IU a. 


EN e o 


2pe sin $ E 


~ a(1—é) 


II, (¢, t)- 7 eG 57. 


_ dip pesing dX, pesing 
saa dt a(1—e) dt * ad —e) a LG )-TS 9, 


and substituting in these equations the values of 


a EIE A and 5 d j^ V(1 — e) 


we find 


T=f2 Poos(v,—A)—1+ c y (cos (v, —A)-— 3 va i de, 2 4-28 sin(v—2) aaa T 


pe sin ġ ; p dt 
TML Gat) a (1— e) IG t) ane V(1 — e) dt 


R=- ÍP cos (v, — a)— 1+ c NC n AC METTE —A) Vd dca s 


(TG ee 1G 0, 


which are Hansen’s values, except that dr in the coefficient of II’ (£, t) has been vre- 


dT 
placed by its value. 


2, Stone Buildings, 31st March, 1855. 
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